Annealing of supersolidity in plastically deformed solid ^He 
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We present a numerical study of a continuum plasticity field coupled to a Ginzburg-Landau model 
for superfluidity. The results suggest that a supersolid fraction may appear as a long-lived transient 
during the time evolution of the plasticity field at higher temperatures where both dislocation climb 
and glide are allowed. Supersolidity, however, vanishes with annealing. As the temperature is 
decreased, dislocation climb is arrested and any residual supersolidity due to incomplete annealing 
remains frozen. Our results provide a resolution of many perplexing issues concerning a variety of 
experiments on solid ^He. 



The first reports of non-classical rotational inertia 
(NCRI) in solid ''He, jU |2| interpreted as evidence for 
a "supersolid" phase, triggered an avalanche of addi- 
tional work [3] which largely confirmed the existence of 
an anomalous NCRI signal in the form of a sharp drop in 
the period of a torsional oscillator (TO) filled with solid 
^He. However, the nature and origin of the TO anoma- 
lies have been very much in dispute, especially after the 
experimental observation [1] of an increase in the shear 
modulus of solid ''He at the onset of the NCRI. The pe- 
riod anomalies can be due to the appearance of a new 
supersolid phase (producing an NCRI fraction, NCRIF, 
which decouples from the TO), to changes in the elas- 
tic properties, or to a combination of both. It is not at 
all easy to disentangle these two effects. At one extreme 
it has been suggested that the phenomena involved are 
largely or exclusively elastic anomalies "F. On the other 
hand, experimental evidence for the mass flow equivalent 
of a fountain effect [B] , which is in principle the gold stan- 
dard for the existence of a superfluid field, indicates the 
presence of superfluidity. Recent experiments [7j on the 
effects of dc rotation on the observed NCRIF also suggest 
the occurrence of superfluidity in solid ^He. 

In our view, the above alternatives may be a false di- 
chotomy. The apparent close connection between the 
elastic properties of the solid and those of the puta- 
tive superfluid field may be understood if one assumes 
that extended crystal defects, such as dislocations and 
grain boundaries, play a crucial role in the observed 
phenomena. A prominent role for crystal defects [51 H] 
is indicated by several experimental and theoretical re- 
sults. The NCRIF reported varies dramatically from 
sample to sample [8], and can decrease substantially on 
annealing[S]. Also, quantum Monte Carlo [TOl E] cal- 
culations have shown that superfiuidity can occur along 
dislocation cores and near grain boundaries. As posited 
in Ref. 1121 the apparent contradictions may be explained 
by assuming that the motion of dislocations, known to 
affect the elastic properties of the solid, also has a strong 
effect on the occurrence of superfiuidity along their cores. 
The reported mass fiow through solid ^He has been 



attributed to quantum superclimb arising from the flow 
of atoms though superfluid dislocation cores [13] . 

In this Letter we present the results of a calculation 
that points to a reasonable resolution of these dilemmas 
by modeling the dynamics of a large number of disloca- 
tions using a recent formulation of continuum plasticity 
theory [T4HT9] coupled to a complex scalar field ip [20l - [22] 
that describes superfiuid order. Our argument is as fol- 
lows: It is known ^16j that an initial smooth distribution 
of dislocations spontaneously coarsens into defect-free re- 
gions interspersed with shock-like structures of high dis- 
location density or "cell walls" . Internal stresses at such 
cell walls, which may be large initially, eventually an- 
neal out producing stress free-grain boundaries (repre- 
senting discontinuities only in crystal orientation) at late 
times [T71 [TH] . If couples only to volumetric stress an 
(the trace of the stress tensor), as in [21 , transient su- 
perfiuidity may exist at cell walls with large an as long 
as dislocation climb is allowed. Once cell walls evolve 
into symmetric grain boundaries, this superfiuidity van- 
ishes, in agreement with quantum Monte Carlo results 
[lOj . If climb is arrested, the plastic current is volume 
preserving |17|,ll8j and tp, if initially absent, cannot form. 
Similarly, annealing of an initially non-zero ip is severely 
constrained without dislocation climb. This suggests the 
following scenario for the occurrence of NCRI: At high 
temperatures, dislocation climb, ensured through mass 
transport through dislocation cores [13^, results in the 
formation of symmetric grain boundaries and vanishing 
of the supersolid fraction at long times. As the temper- 
ature is reduced and/or pinning of dislocations by impu- 
rities becomes effective, climb is suppressed and the re- 
sulting long-lived residual supersolidity at cell walls con- 
tribute to the NCRI. 

The Ginzburg-Landau free-energy for the complex su- 
perfiuid order-parameter V'(r) coupled to the elastic dis- 
placement U; of the solid is [201 121] , 



dr 



do 
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The parameters co,ao and dp a-U are temperature depen- 
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dent and the interaction energy is 



(2) 



where = ^{diUj + djUi) and gij couple the elastic 
degrees of freedom to the magnitude of ip. To the extent 
that NCRI is caused by supersohdity of "^He, the norm 
IV'P of the order parameter is expected to be related to 
the NCRIF measured in experiments. 

As a consequence of the coupling in ([2]), the solid ex- 
periences an "external" stess due to ip, 



(3) 



Dislocations within the solid respond to both af^ and 
to the stress due to the presence of other dislocations. 
To address the problem of defect dynamics in a solid, 
we follow the recent work of Acharya [13] and that of 
Limkumnerd and Sethna |15H19| . It is convenient to in- 
troduce the plastic deformation field [H] ^S^-, such that 
the total deformation gradient djUi = (3fj + is a sum 
of elastic and plastic parts [HI. The dislocation density 
tensor is Py(r) = J2'ytlb]S{r - r'') ^ StinidPl^^/dri, 
where the vectors f and are the tangent and Burg- 
ers vectors respectively of a dislocation line at and 
Eiim is the antisymmetric tensor. Dynamical equations 
for /3^- can now be derived after writing the current 

Plj = E7 Jtj = J2'y^iinit]b]vZ,S{r - rT) as a sum over 
single dislocation contributions, with the velocity of a 
single dislocation line Vi — —Deijktjbiai^i proportional 
to the total stress aij — afj + crfj consisting of separate 
contributions from other dislocations (cfj) and from ip; 

is a material dependent time scale over which plas- 
ticity anneals. Apart from being driven in the direction of 
b by the local stress (glide), dislocations may also climb 
i.e. move in the perpendicular direction in response to 
the local flux of point defects. Total volume is preserved 
by glide, but inclusion of climb removes this constraint. 
The difference between glide and climb motion was in- 
corporated phenomenologically in Ref. by writing the 
total flux as a sum of two terms, 
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7,Sij'^kinit]ble„ipq(Tprt^b'J6{r - r'^)] (4) 



For A = 1, Jij is traceless i.e only glide motion is possible 
and for A = glide and climb are equally probable [l4] . 
Next, coarse graining the current over all dislocations, 
within a mean-field approximation, one obtains 
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(5) 



which is identical to that obtained in [16j except for a 
redefinition of aij . The stress due to dislocations is af^ = 



-Cijfcm/Jfcm where, for an isotropic solid, 

Cijkm — f-^i^ik^jm ^" ^ira^jk ^" ~^ij^km) 



(6) 



and /i and v are the shear modulus and Poisson's ratio 
respectively. In Ref. [T6| it was found that an equation 
similar to ([5| spontaneously leads to dislocation pile-ups 
with associated stress jumps [TB], known as cell walls, at 
finite time. At later times, cell walls evolved to symmet- 
ric stress-free grain boundaries |l7j by attracting more 
and more dislocations to themselves, as in the formation 
of finite time shocks in Burgers turbulence |23j . 

In principle, Eq. ([s]) together with a dynamical equa- 
tion for "0 such as the time dependent Ginzburg-Landau 
(TDGL) equation dtip — ~rSH/6ip, is sufficient to de- 
scribe the dynamics of a supersolid with plastic de- 
formation. The parameter F sets the time scale for 
the evolution of ip. To make explicit calculations, we 
now introduce two simplifications. Firstly, since stress 
aij — dT-L/duij relaxes much faster than the dislocation 
configuration, we may assume that for all times, the di- 
vergence of the total stress vanishes, i.e, djaij = and 
the solid is in mechanical equilibrium. Secondly, we re- 
port calculations only for situations where all the fields, 
/3^-, Tp etc. are functions only of one dimension, z describ- 
ing flat cell walls and grain boundaries [TTj dH] . Using 
af'j and afj in the mechanical equilibrium condition, we 
obtain, in one dimension. 



dz 



1 - V 



Plk- 



(7) 



where (i^^. = I3\-^ + /Jfa- We have simplified Tiint by 
assuming that ip couples only to the volumetric stress 
an = [2/i(l -|- i/)/3(l — 2v)\dzUz. A similar assumption 
was used by Toner [5T] with the sign of g chosen such that 
compressive stresses lead to superfluidity. The Hamilto- 
nian is now given by, 



T~L — Hlg + T~Lint — dz Co 



(8) 



with renormalized parameters as a consequence of Eq. ([t]) 
viz., a — ao—g{i'/l — v)(i^^ and d — d^— g'^[{l — 2v) / ^{l — 
v)]. The constants oq and do are both > 0, so that the 
equilibrium bulk solid does not show superfluidity at any 
temperature. The time evolution equation for /J^^- (i,j = 
1, 2) and (3^^ including superfluidity is given by. 
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FIG. 1. (Color-online) (a) Peach-Kohler force density J^^^ 
(blue/dark gray curve) in the presence of supersolidity for 
A = at reduced time t = 50 as obtained from numerical 
integration of Eqs. ([9| and ( |10[ |. The red/light gray curve 
shows in the absence of ^ (i.e. (7 = 0) for the same initial 
conditions for comparison, (b) plasticity tensor components 
/3jj, /322 from the same calculation as in (a). 
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FIG. 2. (Color onhne) Time evolution for t < 50 of (a) + 
/322 = /^fcfc and (b) for the system with parameters as in 
Fig[l] The supersolid fraction becomes large at one of the cell 
walls where is also large (and so is the absolute value of 
the volumetric stress an). It subsequently decays to zero as 
the stress relaxes. 



with 8 = —<^ijf3fj being the elastic energy density from 
dislocations. For A = 0, /Jgg does not evolve and for 
A = 1, the evolution is constrained such that the trace 
/^ii + 1^22 + P33 vanishes at all times. We choose cq and 
as our units for length and time respectively, the 
energy scale is set by ksT. 

We shall first consider the case with A = 0. The com- 
ponents of /3f (see Eq. (IqI)) are all coupled through the 



Peach-Kohler force density (PKF) J"^^ = d£/dz. It 
is instructive, therefore, to obtain the time evolution of 
. For A = we obtain from Eq. (|9|). 



dtJ" 



PK 



PK 
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ga,(l^pj-^^a,/3P,) = o. (10) 



For 3 = Eq. (10) becomes the Burgers equation j23j. 

We have solved the coupled Eqs. ^ and ( 10 1 together 0.4 
with the TDGL equation for -0, numerically using an ac- 
curate fourth order Runge-Kutta scheme [24j with time 0.3 
step At — 0.05 and spatial discretization Az = 1. The 
initial input for (i^- is random, drawn from a Gaussian 
distribution with mean zero and width — 0.15. The equa- 
tions are regularized by adding a diffusive term ad'^fSf^ °^ 
with a small initial value of a which is subsequently re- 
duced further at later times. The rest of the parameters 
fi — 1, ly = .49, oo — 0.01, g — .5 and d — 1 are cho- 
sen to represent a generic solid above the bulk superfluid 
transition. Finally, for the dynamical parameters we use 
T = D which represents a scenario where \^p\'^ relaxes 
together with the plasticity and the associated stress. 

In Fig. [ija) and (b) we plot, respectively, the PKF 
and the two diagonal components of f3f- . Both the PKF 



and the /3f^ have two discontinuities at z = 180 and « 300 
which are inherited from the Burgers-like terms in Eq. ([9| 
modified by the presence of -0. Since the dislocation den- 
sity pij is given by z derivatives of f3fj, these regions of 
Pfj are also regions of large dislocation content signify- 
ing the presence of dislocation pile-ups or cell walls [IB] . 
While the first of these cell walls is, coincidentally, almost 



free of compressive stresses, the second has a prominent 
stress jump which reduces a in Eqn.([8| locally. The loca- 
tion and nature of the discontinuities depend, of course, 
on the realization of the random initial condition used. 

In Fig. and (b) we show, respectively, the time 
evolution of and of at these cell walls. Initially, 
there are small pockets of superfluidity at random val- 
ues of z, where the local stress is suitably compressive so 
that a < 0. As a consequence of the plasticity dynam- 
ics, dislocations bunch together producing cell walls after 
some time i ~ to- Supersolidity is associated only with 
those cell walls which have a positive stress jump. While 
this condition is satisfied at the second cell wall, the first 
one remains free of ip- Subsequently, the stress decays as 
^ l/\/t — to [m |23] driving the local a > 0, causing ip 
to vanish after some time. 
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FIG. 3. (a) Magnitude of the order parameter at t = 50 for 
A — 1.0, A — 0.6 and A = 0.0. For A = 1, formation of 
regions with large \jp\'^ is suppressed, (b) Spatially averaged 
\tf}\^ vs time t for protocols A and B (see text). A: A = with 
other parameters identical to Fig[2]and Fig[T] B: here A was 
increased to 1 after t = 12.5, note the arrest of the vanishing 
of IV'P in this case. 



If "0 couples only to the volumetric stress as in our case 
(and in Ref. [HJ, then without dislocation climb, we do 
not find supersolid behaviour. This is shown in Fig. [3j^a) 
where we have plotted for different values of A. For 
A = 1, when climb motion is completely suppressed (a 
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reasonable limit for solids at low temperatures far from 
melting) , supersolidity never appears since J^- is traceless 
(see Eq. Q) and compressive stresses are not sponta- 
neously produced. On the other hand, any pre-existing 
compressive stress relaxes extremely slowly due to the 
constraints in Eq. This suggests a strong link be- 

tween suppression of dislocation climb and supersolidity 

In Figjsjjb) we illustrate this link. We show the evolu- 
tion of the spatially averaged jf/'p as a function of time 
following two protocols A and B. Protocol A is identi- 
cal to that shown in Figj2] with A = 0. On the other 
hand in B, we increase A to unity after simulating with 
A = up to < = 12.5, mimicking a decrease in the rate 
of climb. Figjsjb) shows that, in this case, the anneal- 
ing of tp is arrested and cells walls with large may 
persist at low temperatures to experimentally observable 
times. Finally, the experimental time over which NCRIF 
would become unobservable in solid "^He does, of course, 
depends on the parameters F and D. If F > D and the 
plasticity relaxation time, D~^, large, NCRI phenomena 
will be ubiquitious but may depend sensitively to sample 
preparation history [El [9]. Further, D is, in general, also 
dependent on the dislocation density pij (work harden- 
ing) and pinning by ■^Hc impurities [25j at low temper- 
atures, both of which will act to arrest the annealing of 
ip and strengthen our conclusions. However, inclusion of 
such effects may require going beyond the version of con- 
tinuum plasticity used here [16l [19] . On the other hand 
if F << D, NCRI phenomena may become altogether 
unobservable since the required stress relaxes before a 
non- vanishing ip can develop. 

The results obtained from our coarse-grained descrip- 
tion of superfluidity and plasticity are consistent with 
those of earlier microscopic studies [101 US] using quan- 
tum Monte Carlo methods. In Ref. (TUl it was found that 
grain boundaries separating crystallites that differ from 
each other by a simple rotation do not support superflu- 
idity. Other, more generic grain boundaries were found 
to exhibit superfluidity in regions of a few lattice spacings 
in width. Such microscopic regions of superfluidity will 
not appear in our coarse-grained description. Supersolid 
behavior may, of course, arise from a system-spanning 
network of superfluid channels [201 HH HI] of microscopic 
width (e.g. cores of dislocation lines fVD and microscopic 
regions near grain boundaries) . A description of this be- 
havior is beyond the scope of our coarse-grained analy- 
sis. The importance of dislocation climb in supersolidity, 
found in the present work, has been emphasized earlier 
in the microscopic study of Ref. [T31 

Our work shows that supersolidity should be observ- 
able in solid *He when certain thermodynamic, mechan- 
ical and dynamical conditions are satisfied. We find that 
the degree of supersolid behavior depends crucially on 
the details of the annealing process. This necessarily im- 
plies a large sample to sample variation, as observed in 
experiments [8] . Our observation that supersolidity van- 



ishes in the long-time limit if the plasticity field evolves 
according to its natural dynamics (including dislocation 
climb) provides a reasonable explanation of the experi- 
mental observation ^4^ that the shear modulus of solid 
^He increases as T is lowered across the onset tempera- 
ture of the NCRIF. The increase in the shear modulus is 
generally attributed [4, to the pinning of dislocation lines 
at impurities such as "^He atoms [25] . Pinning should 
suppress climb motion, thereby preventing the annealing 
out of the supersolid field. So, the onset of supersolid 
behavior and pinning-induced enhancement of the shear 
modulus should coincide in this picture. Our results on 
the importance of climb motion in the initial development 
of the supersolid order parameter provide an explanation 
of experimental results indicating that mass fiux, dislo- 
cation climb and NCRI are coupled together in solid "'He 
[H HTl [13] . Extensions of this work to investigate the role 
of ^He impurities f5S] and to obtain the elasto-plastic re- 
sponse of a supersolid with dislocations to external stress 
and the glass transition [2 7) in this system are in progress. 
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